Preprint HU-TFT-93-44 



BUBBLE GROWTH AS A DETONATION 

M. LaineB 

Department of Theoretical Physics 
P.O. Box 9, 00014 University of Helsinki, Finland 

7 September 1993 



Abstract 

In the case of spherically symmetric chemical burning, only very special hy- 
drodynamical detonation solutions exist. These are the so called Jouguet det- 
onations, in which the burning front moves at sound velocity with respect to 
the burnt matter. Usually it is believed that the situation is similar in the case 
of bubble growth in cosmological phase transitions. In this paper it is shown 
that actually a much larger class of detonation solutions exists in cosmological 
phase transitions. 
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1 Introduction 



First order phase transitions in cosmology — probably at least the electroweak phase 
transition and the quark-hadron phase transition — proceed through the supercooling 
of the high-temperature phase and the nucleation of supercritical bubbles of the low- 
temperature phase. The nucleated bubbles start to grow. Assuming that this process 
is mainly hydrodynamical, the bubbles can grow either as deflagrations or as detona- 
tions. In a large part of parameter space, the bubbles grow as deflagrations which 
must be preceded by shock waves to satisfy the boundary conditions. However, if the 
supercooling is considerable the bubbles can grow as detonations [§, |fl. 

Both deflagration and detonation solutions can be further subdivided into three 
categories which have different qualitative features [[l], |2|]. The categories are called the 
strong, the weak and the Jouguet detonations and deflagrations. These different types 
of processes are characterized by the velocity of the combustion front with respect 
to the matter behind the front, the Jouguet processes being those for which matter 
flows out of the interface at sound velocity. In any attempt to account for the precise 
mechanism of a first order phase transition one should in the very beginning be aware 
of the category to which the growing bubbles belong. However, to be able to do so one 
must give up the idealization of viewing the phase transition surface as a discontinuity 
and to study its microscopic physics. In the context of chemical burning this problem 
is well understood |^|, In the case of cosmological phase transitions, less is known. 

In this paper we show what the relevant categories of detonations are, and we 
also review the situation in the case of deflagrations. For chemical burning taking place 
as a deflagration, it is well known that strong deflagrations are impossible. This means 
that the burning front cannot move supersonically with respect to the burnt matter. 
This result will be restated in section 4 in the context of cosmological phase transitions. 
However, nothing more can said about the velocity of the phase transition surface with- 
out microscopic calculations such as those performed in refs. for the electroweak 
phase transition and in ref. || for the quark-hadron phase transition. For bubbles 
growing as detonations the situation is different and quite interesting. Namely, in the 
case of chemical burning the bubbles can only grow as Jouguet detonations [[[], || . This 
remarkable result (called the Chapman- Jouguet hypothesis) means in particular that 
the velocity of the combustion front is completely determined by energy-momentum 
conservation and by the boundary conditions, and no degrees of freedom are left to 
the microscopic physics of the burning surface. In the case of phase transitions, there 
has been some confusion about the validity of this result. Steinhardt has argued in 
ref. || , using the same arguments as Landau and Lifshitz used for chemical burning in 
ref. fll], that the Chapman- Jouguet hypothesis should be valid. His proof was originally 
in 1+3 dimensions but it is very easily "extended" to 1+1 dimensions as well. The 
authors in for instance refs. || |] were also under the impression that only the Jouguet 
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detonations are possible, though at least in refs. @, [10| the weak detonations in 1+1 
dimensions are accepted as well. There seems to have prevailed the belief that the 
dimension of the space would have something to do with the validity of Steinhardt's 
results. In refs. IpTTL |12f the authors use Steinhardt's results in the 1+3-dimensional 



context, and also in ref. |13j there is some uncertainty about the possibility of the 
different kinds of detonations. It is the purpose of this paper to show that due to 
qualitative differences between chemical burning and phase transitions, the Chapman- 
Jouguet hypothesis does not hold in the case of phase transitions and therefore weak 
detonations are possible, in addition to the Jouguet detonations. 

The plan of the paper is as follows. In section 2 the relativistic hydrodynamics 
of first order phase transitions is reviewed. In section 3, it is proved that strong 
detonations are not possible as a means of bubble growth in phase transitions. In 
section 4, it is shown that strong deflagrations are impossible as well. In section 5 we 
address the question of weak detonations and show that these are a natural means of 
bubble growth in phase transitions. The conclusions are in section 6. 



2 The hydrodynamics of bubble growth 

Consider a spherically expanding combustion front. Locally the combustion front looks 
planar. If the front velocity is constant, we can Lorentz transform into the rest frame 
of the front. In this frame, we denote the matter ahead of the front (the unburnt or the 
"quark" matter) with subscript 1 and the matter behind the front (the burnt or the 
"hadron" matter) with subscript 2. It is a very good first idealization to treat the front 
as a discontinuity. The strong deviation from thermodynamical equilibrium attached 
with the burning of the quark phase to the hadron phase and all the entropy production 
is confined to this discontinuity, and outside it the energy-momentum tensor is that of 
an ideal fluid: T^ u = wu^u" — pg^ v . Here w = e + p is the enthalpy density, e is the 
energy density, p is the pressure and = (7, 7V ) is the four- velocity of the fluid. The 
energy-momentum conservation across the discontinuity yields the equations 

W\l\v\ + pi = W 2 llv\ + p 2 ■ 

The non-negativity of entropy production is expressed as 

^272^2 > S17JV1 , (2) 

where s is the entropy density. For generality, we also assume that there is a conserved 
quantity (its density is denoted by n), which is in thermodynamical equilibrium with 
the rest of the fluid. Then we have the additional equation 

mjtV! = n 2 ^ 2 V2 = j ■ (3) 
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In the quark-hadron phase transition the conserved quantity is the baryon number. In 
the electroweak phase transition no such conserved quantity exists, but one can always 
add a fictitious conserved quantity (a "tracer" ) to the fluid as long as it does not show- 
in the equation of state. To see how this can be done, consider a problem in which 
there is no conserved quantity, \i = 0. Then the equation of state is p = p(T) and the 
complete solution of the hydrodynamical problem (possibly involving discontinuities) 
consists of the functions T(t, x), v(i, x) and <f)(t, x). Here <j)(t, x) is the order parameter 
and it is needed when the problem contains a phase transition. At the initial time t — 0, 
choose an arbitrary function n(0,x), for instance n(0,x) = 1 in some units. Then one 
can integrate the first order partial differential equation d^nu^ 1 ) = 0, where is now 
a known function of t and x. The solution n(t, x) is the required conserved quantity. 
Although /i = 0, the quantity n enters our thermodynamical equations, because it is 
natural in the present context to formulate the thermodynamical identities using the 
entropy per tracer a = s/n. Then the quantity x = w/n 2 will prove to be very useful, 
see below. Notice also that the arbitrariness of the initial condition for n(t, x) means 
just that it is equivalent to say that the entropy per one tracer is increased, or that the 
entropy per ten tracers is increased. Finally, in a thermal fluid with vanishing chemical 
potential, the true physical particle density n of all the massless relativistic particles 
(and antiparticles) is proportional to entropy and therefore satisfies the inequality 
d^ifiu^) > in contrast to the equation satisfied by the tracer. The inequality sign is 
obeyed at least at the discontinuities. In the following, when we speak of baryons we 
mean just some conserved quantity, possibly the tracer. 

If one does not want to use the tracer in the case \i = 0, it is still possible to 
repeat the following analysis almost as such, see ref. The quantity x only has to 
be defined in another way. 

To describe the state of the fluid on both sides of the discontinuity, two intensive 
thermodynamical variables are needed. It is conventional to take as these the pressure 
p and the variable x = w/n 2 . Then the state of the fluid is represented by a point in 
the (x,p) plane, and the above continuity conditions relate the two points. Explicitly, 
from equations ([[]) and @ we get the equations 

~f = (P2-Pl)/(X2-Xx) (4) 
»2^2 - Will = (p 2 -Pl){x2 + Xi) . (5) 

Given the point (xi,pi) and the equations of state of both phases, the enthalpy w can 
be expressed in terms of x and p and from equation (^J) the curve p 2 = ^2(^2) can be 
solved. This curve is called the detonation adiabat. If both phases are ascribed the 
same equation of state, then eqs. (f|) and (^) do not describe a combustion front but 
a shock front within one phase. In this case the curve p 2 = ^2(^2) is called the shock 
adiabat or the Taub adiabat. 

To illustrate the general structure of the detonation and the shock adiabats in 
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the (x 2 , p2)-plane, we need to know something about the equations of state of the 
burnt and the unburnt matter. In view of the comments after eq. (H), we could use for 
instance the bag equation of stateS in this purpose as was done in ref . || , although this 
would mean that the chemical potential vanishes and hence there is no true conserved 
physical quantity. However, we can easily work a bit more generally to understand some 
of the underlying assumptions. Consider first the shock adiabat of the unburnt matter 



through the point (xi,pi). Taking either the equation of state in ref. [|15|, eq. (2.3)] 
with non-zero chemical potential, or the bag equation of state, we obtain the relation 
w\ = 4pi + 4B. Then from eq. (H) the equation of the shock adiabat becomes 

1 18 
[Ps + 3 (Pi + 45)] [x s - 3 *i] = gZi (Pi + B) . (6) 

The curve p s = p s (x s ) is a hyperbola and it is drawn schematically in figure 1 with 
dashed line. Consider then the detonation adiabat. It is essential that the detonation 
and the shock adiabats never cross and that the detonation adiabat always lies above 
the shock adiabat. To prove the first claim, notice that it follows from equation (|5]) 
that if the two adiabats crossed, then the final states would have exactly the same 
enthalpy density, particle density and pressure. This is not possible because the two 
phases have a different equation of state: w\(n 2 ,p 2 ) ^ W2{n 2 ,p 2 ). For instance, at 
the critical temperature in the quark-hadron phase transition we have s\ (T c , fi c ) > 
s 2 (T c , n c ) because of the latent heat and rti > n 2 as is seen in ref. [|15|, fig. 3], so 
that W\ = T c s\ + ii c n\ > w 2 . Specifically, for the bag equation of state we have 
w\{p 2 )—w 2 {p 2 ) = 4-B. To prove the second claim, consider equation @) in the limit that 
X2 approaches X\ from below. It is seen that p 2 > P\ for all x 2 < X\. Therefore, because 
the adiabats cannot cross even at x 2 — x±, the point (xi,p 2 ) must lie above (xi,Pi) 
and the claim is proved. With this knowledge, the general structure of the detonation 
adiabat can be drawn. In many cases (e.g. in the bag equation of state in the limit 
B — > 0) it is even so that tuning some parameters makes the equations of state of the 
burnt and the unburnt phases the same so that near this limit the detonation adiabat 
resembles greatly the shock adiabat. To prove more precise statements about the 
detonation adiabat, one should combine the non-relativistic analysis in refs. fl], ||| with 



the relativistic analysis in ref. [16|. The detonation adiabat is represented schematically 
in figure 1 with solid line. 

We can now enumerate the different qualitative processes by which the bubble 
growth can proceed. From equation (|j) it is known that the line in the (x,p) -plane 
through points (xi,p\) and (x 2 ,p 2 ) always has a negative slope. Therefore, there are 
two kinds of solutions: either p 2 > p x or p x > p 2 . Solutions of the first kind are 
detonations; examples of these are the points B and C in figure 1 where the dash-dotted 



the bag equation of state, the pressures and the energy densities are pi = a\Ti — B, e\ — 
3aiT* + B , p 2 = CL2T2 and e 2 = 3a 2 T2- Here a\ = a 2 + B/T^ and B is the bag constant. 
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line intersects the detonation adiabat. Solutions of the second kind (points F and G in 
figure 1) are deflagrations. There are three kinds of both detonations and deflagrations. 
Point C is a strong detonation, point B is a weak detonation and point E is a Jouguet 
(or Chapman- Jouguet) detonation. Similarly, point G is a strong deflagration, point 
F a weak deflagration and point H a Jouguet deflagration. From energy-momentum 
conservation alone none of these processes can be excluded. Notice, however, that the 
entropy condition (0) does in many cases rule out the whole family of detonations 
fig. 17] but it does not categorically rule out any of the different types of detonations: 
when the phase transition is preceded by sufficient supercooling, all the processes are 
in principle possible. 

The different types of processes can be characterized by the velocities at which 
matter flows into and out of the discontinuity. To be able to do so, note first that the 
slope of a shock adiabat through any point in the (a;, p) -plane is the quantity —n 2 j 2 v 2 , 



where v s is the local sound velocity ||T6|. Second, consider the detonation adiabat 
through the point (x2,P2), which is hereafter called point 2. This curve is not the same 
as the shock adiabat through point 2. It can be seen using eqs. (HI) and (|5l) and the 



relativistic versions in ref. of the results of ref. []T], §129] that the detonation and the 
shock adiabats through point 2 cross at exactly two points, namely the points B and 
C in figure 1 when point 2 is either of these, and that the slopes of the detonation and 
the shock adiabat agree only if point 2 is point E. Therefore, contrary to Steinhardt's 
claim in ref. ||, the slope of the detonation adiabat does not give the local value of the 
quantity —n 2 j 2 v 2 . However, when the slope of the detonation adiabat is bigger than 
the slope of the straight line between the points (xi,Pi) and (x 2 ,P2), then also the slope 
of the shock adiabat at point 2 is bigger than the slope of this line, and vice versa. 
Because the slope of the straight line between (x\,p\) and (0:2,^2) is by eq. (HJ) just 
— j 2 = —n\y[v\ = — n^fff , we see directly from figure 1 the following characteristics 
of the different processes: 



strong detonation: v\ > v s i, V2 < v S 2 
Jouguet detonation: V\ > v s i, v 2 = v s2 
weak detonation: V\ > v s i, v 2 > v s2 
strong deflagration: V\ < v s i, v 2 > v S 2 
Jouguet deflagration: v\ < v s i, V2 = v S 2 
weak deflagration: V\ < v s i, t> 2 < v s2 



Here v s \ is the local sound speed of the unburnt matter at the thermodynamical state 
(xi,pi) and v S 2 is that of the burnt matter at (22, ^2)- -^ or instance, a weak detonation 
is supersonic relative to the matter both behind and in front of the surface. To relate 
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the velocities v\ and v 2 to the propagation velocity of the combustion front, boundary 
conditions must be considered. The boundary conditions of an expanding bubble are 
that the matter far ahead of the phase transition surface (where no information of the 
phase transition has yet arrived) and far behind the phase transition surface (inside the 
bubble of the low temperature phase) is at rest. To satisfy these boundary conditions, 
the flow profile of a growing bubble consists of several regions [fT7[| . In the case of 
a detonation, the matter ahead of the combustion front is at rest (see e.g. fig. 2). 
Therefore the velocity of the detonation front if<j e t is the velocity v\ at which matter 
flows into the combustion front. In the case of a deflagration, the matter behind the 
combustion front is at rest (see e.g. fig. 3). Therefore the velocity of the deflagration 
front fdef is the velocity v 2 at which matter flows out of the combustion front. Hence, 
for instance, weak deflagrations move subsonically. 

It will be useful to notice that formally, treating every front as a discontinuity, 
a detonation front is equivalent to a shock front which is followed by a deflagration 
front moving with the same velocity as the shock front. We prove this for planar 
fronts. Consider a process in which the state of the fluid changes from (x±,pi) to the 
intermediate state (x ,p ) and then from (x ,p ) to (x 2 ,p 2 ). From equation (|5p we get 

w x - WiZi = (p - Pi)(x + an) 
w 2 x 2 - w x = (p 2 - Po)(x 2 + x ) . 
Summing these together we see that the points (an,Pi) and (x 2 ,p 2 ) satisfy eq. (H) if 

PQ-Pl _ P2 ~Pl 

Xq — an x 2 — an 

Therefore, the point (xq,po) is just the point D in figure 1 and the route AC is equiv- 
alent to the route ADC. The part AD is a shock and the part DC a deflagration. 
From the above it follows that the baryon flux j attains the same value at the shock 
front and at the deflagration front. But then matter flows out of the shock front at the 
same velocity as it flows into the deflagration front, which in the planar case implies 
that the shock and the deflagration fronts move with the same velocity. In a similarity 
solution, which is valid after the early stages of bubble growth, the fronts are therefore 
at the same place. Finally, because matter flows out of a weak detonation front with 
supersonic velocity, a weak detonation is equivalent to a shock and a strong deflagration 
and a strong detonation is equivalent to a shock and a weak deflagration. 

3 The impossibility of strong detonations 

Of the presented six mechanisms of a phase transition, only some are actually realized 
as physical processes. To begin with, an expanding bubble of the low temperature phase 
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cannot grow as a strong detonation. This can be seen even on very general arguments. 
Basically, in a strong detonation there are not enough degrees of freedom to adjust to 
arbitrary boundary conditions 0] and specifically not to those of an expanding bubble. 
We will next prove the impossibility of strong detonations very explicitly, too, since 
this proof will provide us with valuable information. For simplicity, our proof will be 
in 1+1 dimensions. An analogous proof in 1+3 dimensions has been given in ref. |]J 
for the non-relativistic case and with some modifications in ref. || for the relativistic 
case. 

The equations governing the evolution of a bubble are the jump conditions in 
eqs. (^) and (||]), the energy-momentum conservation d^T^ = and the charge con- 
servation d^{nu^) = 0. Denote by a the entropy per baryon. Projecting the equation 
dfjT^ = in the direction of and in the direction perpendicular to the flow velocity 
the energy-momentum and charge conservation equations can be written in the form 



d t (nj) + d x (wyv) = 
d t a + v d x a = 

d t v + vd x v = —(d x p + vd t p)/wj 2 

For a similarity solution depending only on the variable £ = x/t these equations reduce 
to 
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(7) 



lv)p'/wy 2 (9) 

where the prime denotes differentiation by £. From eq. (|7|) it follows that the solution 
cannot be v — £. Then eq. (^) implies that a' = 0. But now using simple thermody- 
namics we can write p' in terms of n': 

dp ( dp\ ( de\ dn 2 . dn 

Then eqs. (0) and (|9|) give an equation for v(£) alone, which is trivially solved to yield 
two solutions, namely v' — and 

„.({) = f=if - . (10) 

Putting these two types of solutions together, we get a picture of the bubble growing 
as a detonation. For < £ < v s , the fluid is at rest. At £ = v s starts the rarefaction 
solution, eq. (|10D . Then there are two possibilities: either the phase transition disconti- 
nuity comes straight after the rarefaction solution, in which case it follows from eq. (|10|) 
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that matter flows out of the discontinuity at exactly the sound speed, or there is an 
area of constant velocity between the rarefaction and the discontinuity. In the latter 
case matter flows out of the discontinuity at a velocity greater than v s because the right 
hand side of eq. (|10D is a growing function of £ for constant v. The former possibility is 
a Jouguet detonation, the latter a weak detonation. These are shown schematically in 
figure 2, together with the corresponding 1+3 -dimensional detonation bubbles. Hence, 
strong detonations are not possible. 



4 The impossibility of strong deflagrations 

Another type of solutions which are not realized in nature is strong deflagrations. 
Performing some causal analysis near the interface, one can conclude that strong de- 
flagrations would not be mechanically stable and hence they cannot exist [p], §131]. 
Another way to exclude these is to give up the idealization of viewing the phase tran- 
sition surface as a discontinuity and to study its microscopic structure. The essential 
feature which makes strong deflagrations impossible is that thermodynamical variables 
- to the extent that they can be used in this region — change monotonically and 
continuously between the two phases 0. Then in a strong deflagration the state of 
matter would slide from point A in figure 1 to point G. Assume that inside the burning 
zone matter can be described as a changing mixture of the low-temperature phase and 
the high-temperature phase, so that the energy-momentum tensor is that of an ideal 
fluid!. Then equation (H) is satisfied at every intermediate point (x 2 ,p 2 ) between A 
and G. This means that at point F the state of the system characterized by x 2 , p2, w 2 
(from eq. @) and h 2 (from the equation x = w/n 2 ) has not only the same variables 
x 2 and p 2 as the low-temperature phase at this point but also the same w 2 and n 2 , 
since the detonation adiabat is just the solution of eq. (|). Because the equation of 
state fixes the relation w 2 = w 2 (n 2 ,p 2 ), the system actually already has to be in the 
low-temperature phase at point F. But then the transition from point F to point G 
would be completely equivalent to a rarefaction shock, since it was noted earlier that 
the shock adiabat through point F crosses the detonation adiabat also at point G. Since 
entropy increase does not allow rarefaction shocks, there cannot be any microscopic 
mechanism by which to move from F to G, ie strong deflagrations are not possible. 
Notice, however, that entropy increase does allow a similarity rarefaction solution, 
which is not a discontinuity. Therefore it seems possible that if a strong deflagration 
could somehow momentarily be forced to exist, it would tend to split into a similarity 
rarefaction solution and a weak deflagration. 

To give an example of a typical deflagration solution, we cite results from ref. [I8| . 



There a simple model is given for first order phase transitions taking place in relativistic 



2 When there is a microscopic order parameter field, one has to make the proof somewhat differently. 
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matter. The model includes one phenomenological parameter T of dimension GeV~ 
which is a sort of dissipation constant and fixes the entropy production at the phase 
transition surface. The dissipation constant is roughly proportional to the collision 
time r c . Small T means large friction and large T small friction at the interface. 
In this model, it turns out that in deflagration solutions the temperature and other 
thermodynamical variables really behave monotonically at the phase transition surface, 
see fig. 3. Accordingly, only weak deflagrations (and as a limiting case of these, Jouguet 
deflagrations) have been found. 

5 Weak detonations 

We now turn to weak detonations. It has been argued in refs. Q for the non- 
relativistic case of chemical burning and in ref. || for the case of phase transitions in 
relativistic matter that weak detonations are impossible. Then the Jouguet detonation 
would be left as the only mechanism for bubble growth as a detonation. This is called 
the Chapman- Jouguet hypothesis. To see how this arises, consider chemical burning. 
The rate of chemical burning increases rapidly with increasing temperature, often as 
exp(—U/T) where U is some constant factor. Therefore, heat must be supplied to 
the unburnt matter before the burning can begin. In deflagrations, this happens by 
thermal conduction from the burning zone where heat is liberated on account of the 
exothermic character of the reaction. This mechanism is, however, so slow that the 
velocity of the burning front is subsonic [0] . To make a detonation, an entirely different 
means of raising the temperature is needed, and this can be achieved by a strong shock 
wave which ignites the burning. The shock front compresses the matter so that the 
temperature and the pressure rise dramatically, and then the burning zone follows in 
which pressure again decreases (the temperature increases even in the burning zone, 
which is possible because the chemical potential is non- vanishing) . This means that 
the formal equivalence of a detonation to a shock and a deflagration, which was proved 
above, would also describe the true microscopic structure of the detonation front. In 
figure 1, this means that matter first jumps from point A to point D and then slides 
from D to C as thermodynamical quantities change continuously in the burning zone. 
Since a shock in which pressure decreases is forbidden by the condition of entropy 
increase, matter remains in C and cannot continue to B. Another way to put this is 
that a weak detonation is equivalent to a shock and a strong deflagration, and as strong 
deflagrations are impossible, so are weak detonations. From these arguments, the 
would-be mechanism of a detonation in chemical burning must be a strong detonation, 
point C. But because even strong detonations were above proved to be impossible, 
the detonation must correspond to point E in figure 1. And this is just the Jouguet 
detonation, in accordance with the Chapman- Jouguet hypothesis. 

The question is now what the microscopic structure of a detonation front in 
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the case of a phase transition actually is. First of all, no preheating such as that 
caused by the shock front in the case of chemical burning is needed to start the phase 
transition. Namely, in the case of phase transitions, there is no local "burning rate" 
which would have to be high enough for the transition to start. The only relevant 
rate is the rate of nucleation. This is of the form exp(—U/T) where U is not constant 
but includes singular temperature dependence: U oc 1/(1 — T/T c ) 2 (in the limit of 
small supercooling ||, [OJ). The lower the temperature is, the faster is the rate of 
nucleation. Once the nucleation has happened, the growth of the bubble is determined 
by the equation of motion of the order parameter and by hydrodynamics. It seems 
highly unlikely that the solution of these equations would include a sharp narrow peak 
in pressure, because due to the vanishingly small chemical potential in the case of 
cosmological phase transitions, the peak in pressure would also cause a sharp narrow 
peak in temperature. Such a peak in temperature is by itself rather unnatural, and 
it would also make the two minima of the effective potential more degenerate and the 
potential barrier between the minima higher. Both factors tend to make the tunneling 
between the minima more difficult. Notice also that in the case of deflagrations, a 
rise in temperature after the shock front is necessitated by boundary conditions and 
hydrodynamics, but the supposed temperature peak at the detonation front serves no 
such purpose. For these reasons, it seems much more natural that the order parameter 
- if one exists — and the thermodynamical variables change monotonically to their 
new values in the transition front and the latent heat released in this process fuels the 
rapid expansion of the detonation front. This means that in figure 1 the state of matter 
changes directly from point A to point B. As an example, we quote results from the 
above-mentioned model of ref. In figure 4, a bubble growing as a detonation is 



shown. At the detonation front, the temperature really changes monotonically and the 
solution, accordingly, is a weak detonation. All in all, weak detonations seem to be a 
natural mechanism for bubble growth in phase transitions. 

To be concrete, let us investigate the quark-hadron phase transition in the early 
universe. The physical parameters describing the first order phase transition are the 
critical temperature T c , the latent heat L, the surface tension a and the correlation 
length l c . We take the values a = 0.1 and l c = 6 in appropriate powers of T c . For the la- 
tent heat we assume a small value L = 0.1T C 4 , which results in considerable supercooling 
so that the nucleation temperature is T f = 0.891T C [|J. With this much supercooling, 
detonations are possible. With the model of ref. |18| and the phenomenological param- 
eter r describing the physics of the phase transition surface, figure 5 results. For small 
T the bubbles grow subsonically as weak deflagrations; for large T they grow super- 
sonically as weak detonations. Let us mention in passing that this set of parameters 
yields a very interesting scenario for the quark-hadron phase transition. The average 
distance between the nucleated bubbles would be very large, the bubbles would grow as 
detonations and the phase transition would leave behind it large-scale inhomogeneities 
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6 Conclusions 

When entropy increase allows it, a bubble of the low-temperature phase can grow either 
as a weak deflagration or a weak detonation (the Jouguet processes are here considered 
to be limiting cases of these). For weak detonations to be possible, the nucleation 
must be preceded by considerable supercooling. In weak deflagrations, the velocity of 
the phase transition surface is smaller than or equal to sound velocity, and in weak 
detonations, it is larger than sound velocity. In neither case can the energy-momentum 
conservation and the boundary conditions alone determine the velocity of the phase 
transition surface: its microscopic physics, for instance the entropy production, must 
be known before the exact expansion velocity can be calculated. This means that the 
detonation solutions are qualitatively different in the cases of chemical burning and 
phase transitions. 
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Figure captions 



Figure 1: A schematic picture of the detonation and the shock adiabats. Point A 
is the initial state (xi,pi) of the unburnt fluid. The detonation adiabat 
shows the possible final states in the (x2,p2)-plane. 

Figure 2: On the left is a schematic illustration of the velocity profiles in a 1+1 
-dimensional weak (solid line) and Jouguet (dashed line) detonation. On 
the right are the same profiles for the 1+3 -dimensional case (17 . 



Figure 3: A deflagration solution at time t = 1800 after the nucleation in the 
1+1 -dimensional model of ref. [I8|. The quantities are measured in 



appropriate powers of T c to make them dimensionless. The front on the 
right is the shock front which has at this time not yet sharpened to an 
exact discontinuity, and the front on the left is the deflagration front 
where the phase transition takes place. Both fronts are moving to the 
right. The velocity of the deflagration front is t>d e f — 0.46. 

Figure 4: A detonation solution at time t = 1800 after the nucleation in the 1+1 



-dimensional model of ref. |lq| . The quantities are measured in appro- 
priate powers of T c to make them dimensionless. The front on the right 
is the detonation front where the phase transition takes place, and the 
front on the left is the rarefaction wave. Both fronts are moving to the 
right. The velocity of the detonation front is t>d e t = 0.72. Although the 
rarefaction wave is very thin in this picture due to the scale of the y-axis, 
it cannot be approximated as a discontinuity, see fig. 2. 

Figure 5: The process by which nucleated bubbles grow in the 1+1 -dimensional 



model of ref. |L8] as a function of the parameter V (in units of 1/T C ) [18] . 
At approximately T = 10 the solution changes from a weak deflagration 
to a weak detonation. 
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